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ABSTRACT 

Few high-redshift, radio-loud quasars are known to date. The extremely luminous, 
radio-bright quasar, SDSS J013127.34—032100.1 was recently discovered at a redshift 
of z = 5.18. We observed the source with high resolution very long baseline inter¬ 
ferometry (VLBI) at 1.7 GHz with the European VLBI Network (EVN) and found a 
single compact radio component. We estimated a lower limit to the brightness tem¬ 
perature of the detected radio component, Tb ~ 10^^ K. Additionaly, when compared 
to archival radio data, the source showed significant flux density variation. These two 
hndings are indicative of the blazar nature of the source. 
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1 INTRODUCTION 

High-redshift {z > 4.5) quasars are of major importance 
since they can provide information on the growth of the 
supermassive black holes and the evolution of active galac¬ 
tic nuclei (AGN) in the early Universe. The mere existence 
of black holes with a few million solar masses (or more) at 
z ^ 6 co nstrains the blac k hole growth and the accretion pro¬ 
cess (e.g. lVolonterillioi^ . and references therein). Among the 
high-redshift quasars, the radio-loud ones constitute a very 
attractive subsample, since their radio jets can be studied 
with the highest angular resolution, via very long baseline 
interferometry (VLBI) technique. 

Typically, radio-loud AGN contain a flat-spectrum core 
and a steep-spectrum jet. However, the higher redshift is, 
the fainter the AGN radio jet will appear if observed at 
a fixed wavelength. Thus, if core-jet AGN constitute the 
same population of objects throughout the redshift space, 
the appar ent “prominen ce" of jets at higher redshifts must 
decrease (lGurvitslll999l ): well pronounced jets at high red- 
shifts must appear less frequent than at low redshifts. 
Additionaly, recent VLBI studies suggest that the naive 
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expectation that relativistically beamed sources (blazars) 
should dominate the high-redshift radio-loud quasar pop¬ 
ulation might not be true. There seems to exist a popu¬ 
lation of steep-spectrum high-redshift (z > 4 . 5) rad io-loud 
quasars (iMoniiian et ahlboo^ : l^ev et al.ll2008l . [201oh . which 
may be very young radio quasars, similar to the gigahertz- 
peaked spectrum (GPS ) sources observed in the local Uni¬ 
verse llFrev et al.l I2OIOI) . On the other hand, a few high- 
redshift, beamed blazars are known as well. Thr ee are iden- 
tifled at z > 5: SPS S J114657.79-1-403708.6 llFrev et al.l 
I2OIOI: iGhisellmi et a.1.1 |2014|). SPSS J1 02623.61-b254259.5 
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2013a: 


Frev et ahl gOlSh . and Q0906-I-6930 


I Romani l2006l : ISbarrato et al. 2013lj) . 

Recently lYi et ahl |(20l3) reported the discov¬ 
ery of a high-redshift, radio-bright quasar, SDSS 
J013127.34-032100.1 (J0131-0321, hereafter). The 

source was first selec ted as a candidate high-redshift 
quasar (IWu et al.ll2^ using the optical-infrared selection 
criteria based on the Sloan Digital Sky Survey (SDSS, 
lEisenstein et al] 

Explorer ( WISE, 
cal spectrum was first measured wit h the Yunnan Fai nter 
Object Spectrograph and Camera llZhang et al.1 l2014l ) in 
2013 November. Later higher resolution optical and near- 


2011^ and Wide- Field Infrared Survey 


Wright et al I I 2 OIOI) photometry. Its opti- 
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infrared spectra were obtained by the Magellan Echelette 
and Folded Port Infrared Echelette (FIRE) spectrographs 
in 2014 Jannary. According to these measurements , the 
source has a redshift of a = 5.18 ± 0.01 jYi et al.||2014|) . 

Using the optical spectra and the empirical scaling re¬ 
lation between the Mg II line width and the black hole 
mass (|McLxire fc DunlotJ l2004l l , lYi et al.l (|2014l 'l estimated 
the total bolometric luminosity and the black hole mass 
of J0131-0321 to be Lboi = (1.1 ± 0.2) x 10"^^ W and 
Mbh = 2.7tn'^ X lO^Mpi, respe c tively . Using the latest rela¬ 
tion of iTrakhtenbrot fc Netz^ ll2012ll . the black hole mass 
can be even higher, 4 x 10 ®Mq. 

According to the Faint I mages of the Radio Sky at 
Twenty Centimeters (FIRST, iBecker et al.l Il995h survey, 
J0131—0321 is radio-bright with a flux density of 33.7 mjy 
at 1.4 GHz. Here we report on the results of our high res¬ 
olution radio interferometric observation of the source with 
the European VLBI Network (EVN). 

Throughout this paper, we use flat A cold dark matter 
cosmological model with Ho = 70kms“^ Mpc“^, = 0.3, 
Ha = 0.7 (the same as used in the discovery paper. lYi et al.l 
l2014h . 


2 VLBI OBSERVATION AND DATA 

REDUCTION 

The exploratory EVN observation of J0131—0321 took 
place on 2014 Dec 2 at 1.7GHz (project code: RSG06). 
Six antennas participated in this e-VLBI experiment: Ef- 
felsberg (Germany), Hartebeesthoek (South Africa), Jo- 
drell Bank Mk H (United Kingdom), Onsala (Sweden), 
Toruh (Poland), She shan (China). In an e-VLBI experi¬ 
ment (ISzomorrJ l2008h . the signals received at the remote 
radio telescopes are transmitted over optical fiber networks 
directly to the central data processor for real-time cor¬ 
relation. This was performed at the EVN software cor¬ 
relator (SFXC, iKeimpema et all l2015h in the Joint Insti¬ 
tute for VLBI in Europe (JIVE), Dwingeloo, the Nether¬ 
lands, with 2 s integration time. The observation was car - 
ried out in phase-reference mode (jBeaslev fc Gonwavlll995l l. 
thus precise relative positional information could be ob¬ 
tained. The phase-reference calibrator, J0123—0348 is sep¬ 
arated by ~ 2° from the target in the sky. Its coordinates 
are right ascension ao = 01^23“'35!77473 and declination 
Jo = —03°48'39”3162. The positional uncertainties are 0.3 
and 0.63 milliarcseconds (mas) in right ascension and dec¬ 
lination, respectivelj0. The observation lasted for 2 h, the 
total time spent on J0131—0321 was nearly 90 min. Eight in¬ 
termediate frequency channels (IFs) were used in both right 
and left circular polarizations. Each IF was divided into 64 
spectral channels. The total bandwidth was 128 MHz. 

T he NRAO A stronomical Image Processin g System 
(AlPS ■ [Greisenlf2003h was used for data calibration IIDiamondl 
Il995h following standard procedures. The interferomet¬ 
ric visibility amplitudes were calibrated using the sys¬ 
tem temperatures and antenna gain curves measured at 
the telescopes. After fringe-fitting performed for the phase 

^ Data are from http://astrogeo.org maintained by L. Petrov, 
rfc_2014d solution 
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Figure 1. 1.7-GHz EVN image of J0131—0321. The peak is 
60.1 mJy beam“^, the restoring beam is 5.7 mas x 2.7mas at a po¬ 
sition angle of 82° and shown in the lower left corner of the image. 
The lowest contours are at ±3(7 noise level (±0.25 mJy bearn”^), 
further positive contour levels increase by a factor of two. 

calibrator, its visibilit y data were exported to difmap 
(jShepherd et al.l Il994l ) for imaging. The standard hy¬ 
brid mapp ing procedure was used, with several cycles of 
CLEANing llHogbomI Il979l l and phase self-calibration, then 
finally amplitude self-calibration. Antenna gain correction 
factors were determined in difmap and applied to the data 
in AlPS. The correction factors on average were below 10 per 
cent. Fringe-fitting was repeated for the calibrator source 
but now taking into account its CLEAN component model 
(representative of its brightness distribution) to reduce the 
small phase variations due to its structure. The obtained so¬ 
lutions were interpolated and applied to the target source, 
J0131—0321. In order to correct for the bandpass, four-four 
channels were discarded at the edges of each IF. Then the 
calibrated and phase-referenced data of the target source 
were imaged in difmap. No amplitude self-calibration was 
attempted for the target source. Since J0131—0321 turned 
out to be adequately bright (see Sect. |3]), we were able 
to fringe-fit the source itself in AlPS. The images obtained 
with and without using the phase-calibrator source agree 
within the uncertainties. In Fig.[T] we show the image made 
from the fringe-fitted visibility data. The weights of the data 
points were set inversely proportional to the amplitude er¬ 
rors (natural weighting). 


3 RESULTS OF THE EVN OBSERVATION 

Our EVN observation revealed a single, compact radio com¬ 
ponent (Fig.[T]). Using the phase-referenced data, we derived 
the position of the brightness peak with the maxfit verb in 
AlPS: right ascension at = 01*'31™27!3473 and declination 
Jt = —03°21'0(y.'0791. We estimate that each coordinate is 
accurate within 1 mas. The dominant contribution to the er- 
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Figure 2. Calibrated visibility amplitudes of J0131—0321 as a 
function of the projected interferometer baseline length. The vis¬ 
ibility points are binned for illustration purposes only, model fit¬ 
ting was performed by using all visibility points. The binned am¬ 
plitudes are weighted averages of the visibility points (see text for 
details), the horizontal error bars represent the width of each bin, 
the vertical error bars are the square root of the weighted sample 
variance. Solid line shows the best fit circular Gaussian model, 
dashed line shows the visibility amplitude of an unresolved point 
source with the same flux density as the Gaussian component. 


ror is the positional uncertainty of the phase-reference cali¬ 
brator. Additional errors come from the thermal noise of the 
interferometer phases, and from the phase-referencing over 
2° in the sky. The obtained position is ~ 117 mas away from 
the coordinates given in the FIRST catalogue. 

To quantify the source brightness distribution, we fitted 
the self-calibrated visibilities with a single circular Gaussian 
model component in difmap. The integrated flux density is 
(64.4 ± 0.3) mjy, the full width at half maximum (FWHM) 
size is (0.79±0.01) mas. T he statistical erro rs were calculated 
following the formulae of lFomalond lll999ll . The angular size 
corresponds to a projected linear size of ~ 4.9 pc in the 
rest frame of the source. The fitted size of the Gaussian 
component is close to, but slightly large r than the minimum 
resolvable size (see lKovalev et aDl2005l ). 

The calibrated visibility amplitudes as the function of 
the projected baseline length are shown in Fig. [2] together 
with the fitted Gaussian brightness distribution model (solid 
line) and a point source model with the same flux density 
as the Gaussian model (dashed line). In the figure, we bin 
the data in {u,v) radius, by putting 4110 visibility points 
in one bin and display weighted averages of the visibility 
amplitudes for illustration purposes only; when fitting the 
model, all visibility points were used. Because of the large 
uncertainty of the visibility measurements on the longest 
baselines, we use the size estimate as an upper limit and 
thus we can derive the minimum rest frame brightness tem¬ 
perature of J0131—0321 as 

Tb ^ 1.22 X 10^"^(l-b2)K, (1) 

where S is the flux density measured in Jy, it is FWHM 
size of the component in mas, z is the redshift, and i/ is 


the observing frequency in GHz. The resulting lower limit 
to the brightness temperature is Tb ^ 2.8 x 10^*^ K. This 
brightness temperature is an order of magnitude higher 
than those of th e other high-reds hift (z > 4.5), non-blazar 
radio-loud AGN dFrev et al.ll2010l), and it also ex ceeds that 
of SPSS J114657.79-1- 403708.6 jFrev et al.ll2010l) . a blazar 
dGhisellini et al.ll2014l ) at 2 = 5.0 by two orders of magni¬ 
tude. The 1.7 GHz monochromatic radio luminosity of the 
source is 3 x lO^^WHz”*^. This is comparable to those of 
th e five 4.5 < z < 5 radio-loud quasars studied by VLBI 
by Frev et al.l ll2 010h and th e three 2 > 5 blazars dl omanil 
l2006l: iFrev et al.ll2010l. l2015l) FI 

Assuming an intrinsic brightness temperature for the 
source (7 b“‘) and knowing a lower limit to its apparent 
brightness temperature, a lower limit to the Doppler fac¬ 
tor (5 = Tb/Tb"*) can be estimated. Often the equip artition 
brightness temperature (5 x 10^° K. iReadheadlllOM ) is used 
as the intrinsic brightness temperature of the source. This 
assumed value indicates moderate Doppler boosting factor 
as a lower limit, 5 > 6. 

The flux density of the source measured in our 1.7-GHz 
VLBI experiment is significantly larger than the value given 
in the 1.4-GHz FIRST cata logue, 33.7mJy and the value 
given in the NVSS catalogue (ICondon et al.lll998h . 31.4 mJy. 
This can be most straightforwardly explained with variabil¬ 
ity: J0131—0321 was almost twice as bright at 2014.9 as 
it was at 2009.2 (the mean epoch of the FIRST observa¬ 
tion) and at 1993.9 (the mean epoch of the NVSS observa¬ 
tion). Alternatively, an extremely inverted spectrum would 
be needed to explain the different flux density values. The 
observing frequency of the FIRST survey, 1.4 GHz, and of 
our EVN observation, 1.7 GHz correspond to ~ 8.7 GHz and 
~ 10.3 GHz in the rest frame of the source, respectively. 
The spectral index {a defined as S' ~ i/“) between these 
frequencies would be ~ 3.8, exceeding the canonical value 
of 2.5 f or a single, optically thick, s ynchrotron self-absorbed 
source (iRvbicki fc LightmanI [l98m . (Spectral index values 
lower than 2.5 straightforwardly explained as the result of in¬ 
homogeneity in the sources.) Additionally, there is no sign of 
such inverted spectrum across the eight IFs of our EVN ob¬ 
servation. The flux density values measured at the different 
IFs agree within the errors. Thus, the more likely scenario 
is that J0131—0321 brightened considerably with respect to 
its FIRST flux density. Such variability is expected from a 
beamed, blazar-type source. The flux density of the phase- 
reference caibrator in our EVN observation, (97.6±0.3) mJy, 
is lower than the flux density values given for the source in 
the NVSS and FIRST catalogues, 148.4 mJy and 125.7 mJy, 
respectively. Thus, it is unlikely that systematic (e.g. instru¬ 
mental) amplitude calibration errors have led to the high 
flux density value of the target source, J0131—0321. 


4 DISCUSSION 

iGhisellini et al.l ll2015l) conducted Swift satellite observation 
of J0131—0321 between 2014 Oct 23 and 2014 Dec 9. They 


^ In the case of Q0906+6930, we calculated the radio lumi- 
nosity using data from the NRAO VLA Sky Survey (NVSS, 
iGondon et al.lll998h catalogue. 
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described the optical-UV-X-ray spectral energy distribu¬ 
tion (SED) of the sonrce wit h a model containing a jet, a n 
accretion disc, and a torus l|Ghisellini fc Tavecchid boOflh . 
They concluded that the viewing angle of the jet is very 
close to the line of sight, 0 ~ 3° — 5°. Our 1.7-GHz, ex¬ 
ploratory e-EVN observation revealed that the source is 
compact with relatively high brightness temperature and 
moderate Doppler boosting factor. Assuming there is no sig¬ 
nificant misalignment between the scales probed via SED 
fitting and the radio-emitting jet observed with VLBI at 

1.7 GHz, using the viewing angle and the lower limit to the 
Doppler factor, we can estimate the Lorentz factor of the jet 

(7): 

7(1-/3cos0)’ 

where (3 is the bulk speed of the jet material given in the 
nnits of the speed of light, c. The res ulting Lorentz factor is 
■ y ~ 3 (and /3 ~ 0.95). In the model of iGhisellini &: Tavecchij 
(I2009lj . the Lorentz factor is tied to the viewing angle (7 = 
d~^), and thus is more than fonr times larger than ours. The 
size of the radio emitting region 4 .9 pc) is less than eight 
times the dissipation radius given bv iGhisellini et ahl (l2015l l 
2 X 10^® m, corresponding to a linear size of ~ 0.65 pc). 
Thus, to explain the existence of the two different Lorentz 
factors, the jet would have to decelerate significantly within 
this relatively small distance. 

Alternatively, since our size estimate is an upper limit, 
and thus the Doppler factor estimate is a lower limit, using 
the Lorentz factor (7 = 13) and th e view ing angle (3° — 5°) 
from the model of iGhisellini et ^ (l2015fl . the Doppler fac¬ 
tor can be calculated to be in the range between 11 and 18. 
Assuming the equipartition brightness temperature for the 
intrinsic brightness temperature of the source, this would 
imply a measured brightness temperature in the range of 
(5.5 — 9.0) X 10*^^ K. For the same flux density that we mea¬ 
sured, this could be achieved if the radio emitting source 
size would be in the range between 0.25 mas—0.4 mas corre¬ 
sponding to 1.5 pc—2.5 pc in the frame of the source. (Or for 
the same source size, the flux density should be ~ 2 —3 times 
larger, which is however in disagreement with our measure¬ 
ment.) 

Finally, we note that the X-ray and radio observations 
were performed relatively close to each other but they were 
not strictly simultaneous, thus flux density variability of the 
source might hinder the comparison of the models deduced 
from the observations. In radio regime, the source variabil¬ 
ity is indeed detected on a longer timescale of a few years 
between our EVN observation and the FIRST data. Apart 
from the SED fitting, the blazar classification of the source 
is also supported by its flux density variability. 

The three previously known z > 5 blazars h ave black 
hole m asses of (2 — 5) x 10® Mq. According to lYi et ahl 
ll2014h . the mass of the black hole in J0131—0321 is ~ 

2.7 X 10 ®M( 7 ), but ca n be as large as ~ 4 x 10®Mp. How¬ 
ever. I Ghisellini et al.1 ll2015h assuming a lShakura fc SunvaevI 
ili97ar disk model, determined an even larger black hole 
mass of (1.1 ± 0.2) x lO^^M©. This larger mass value no 
longer requires super-Eddington accretion to produce the 
bolometric luminosity of Lboi = (1.1±0.2) x lO"^*^ W, the Ed¬ 
dington ratio would be ~ 0.8. The existence of supermassive 
black holes with few billion solar masses at early cosmolog¬ 


ical epochs (e.g. at redshift 2 > 6, iFan et ^l2003ll poses a 
challen ge to Eddington- l imited black hole growth formation 
model. I Ghisellini et al.1 (I 2 OI 5 I I showed that a highly spin¬ 
ning black hole with low accretion disk efficiency is able 
to explain the existence of a black hole with 10^® Mq at 
2 = 5.18. On the other hand, it could also be that the seeds 
of such supermassive black holes had initially been massive 

1 100 — 10®M(7)) and/o r grew via super-Eddington accretion 
Wvithe fc L^l2012lL 


5 SUMMARY 

We observed the newly discovered 2 = 5.18 radio-bright 
quasar, J0131—0321 with EVN at 1.7 GHz. It was detected 
as a single compact radio source. We calculated a lower limit 
to the brightness temperature, Tb ^ 2.8 x 10^^ K. Assum¬ 
ing the equipartition brightness temperature as the intrinsic 
one, the observed value indicates moderate Doppler boost¬ 
ing, 5^6. The flux density of the source in our EVN ob¬ 
servation is almost twice as high as the the value given in 
the FIRST catalogue, measured ~ 5 years earlier. Such flux 
density variability and the relatively high brightness temper¬ 
ature are indicative of the blazar n ature of the sour c e. Thi s 
is in agreement with the results of iGhisellini et ahl ll2015ll . 
who analysed the optical, UV, and X-ray measurements of 
the source. 

VLBI observations at higher radio frequencies (e.g. 
5 GHz) can provide a more robust limit on the source 
compactness and its brightness temperature. Additional 
constraint on the brightness temperature, thus on the 
Doppler factor and the blazar nature of the source, 
can be obtained fro m variability brightness temperature 
llHovatta et al.1120091) provided by radio flux density mon¬ 
itoring of the source. 
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